Nuclear Magnetism and Electronic Order in 13 C Nanotubes 
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Single wall carbon nanotubes grown entirely from 13 C form an ideal system to study the effect 
of electron interaction on nuclear magnetism in one dimension. If the electrons are in the metallic, 
Luttinger liquid regime, we show that even a very weak hyperfine coupling to the 13 C nuclear spins 
has a striking effect: The system is driven into an ordered phase, which combines electron and 
nuclear degrees of freedom, and which persists up into the millikelvin range. In this phase the 
conductance is reduced by a universal factor of 2, allowing for detection by standard transport 
experiments. 



The physics of conduction electrons interacting with lo- 
calized magnetic moments is central for numerous fields 
in condensed matter such as nuclear magnetism 111, heavy 
fermions [|J , or ferromagnetic semiconductors 0,0, [1, [fif . 
Nuclear spins embedded in metals offer an ideal platform 
to study the interplay between strong electron correla- 
tions and magnetism of localized moments in the RKKY 
regime. In two dimensions the magnetic properties of 
the localized moments @, H| depend indeed crucially on 
electron-electron interactions EH EH, EH, EH • In one- 
dimensional (ID) systems such as single wall carbon nan- 
otubes (SWNTs) electron correlations are even more im- 
portant. For metallic (armchair) SWNT they lead to 
Luttinger liquid physics [lj, [l5|, [l6| . Recently, SWNTs 
made of 13 C, forming a nuclear spin lattice, have become 
experimentally available [13, El, EH HH] ■ Motivated by 
this we study here nuclear magnetism in metallic 13 C 
SWNTs. We show that even a weak hyperfine interac- 
tion can lead to a helical magnetic order of the nuclear 
spins (see Fig. [1} coexisting with an electron density or- 
der that combines charge and spin degrees of freedom. 
The ordered phases stabilize each other, and the crit- 
ical temperature undergoes a dramatic renormalization 
up into the millikelvin range due to electron-electron in- 
teractions. In this new phase the electron spin suscep- 
tibility becomes anisotropic and the conductance of the 
SWNT drops by a universal factor of 2. 

The drastic restructuring of the electron wave func- 
tions through the renormalization is very different from 
the case of two @, [1] or three dimensions [l[ where it 
is, in comparison, weak. The same renormalization leads 
to considerable anisotropy in the electron system: The 
nuclear magnetic field spontaneously breaks the spin ro- 
tational symmetry; it rotates in a plane, which we can 
associate with the spin (x, y) directions (see Fig. [TJ. This 
plane is singled out as an easy-plane through the stabi- 
lization of the electron density wave, and electron cor- 
relation functions become anisotropic between the spin 
{x, y) plane and the spin z direction. We illustrate this 
behavior below through the calculation of the electron 
spin susceptibilities. We emphasize that this anisotropy 



is a crucial feature of the SWNT system studied here 
and appears spontaneously due to strong renormaliza- 
tion of the RKKY interactions. This distinguishes our 
system, in particular, from models with built-in easy-axis 
anisotropy [2l| . 

Model. — We assume that the electrons are confined 
in a single mode tp± in the directions perpendicular to 
the tube axis. The nuclear spins / = 1/2 of the 13 C ions 
on a circular cross section have identical overlaps with 
this transverse mode, and so identical couplings to the 
electrons. Through their indirect RKKY interaction over 
the electron gas they are therefore locked in a ferromag- 
netic alignment (see Fig. []}. This RKKY interaction, 
described below, overrules furthermore the direct dipolar 
interaction between the nuclear spins. The latter is very 
small 22], ~ 10 -11 eV, and shall be neglected hence- 
forth. This allows us to treat the nuclear spins as a ID 
chain of large I = INj_ spins, composed of the sum of the 
N± ~ 50 spins around a circular cross section. Due to 
this, Kondo physics, which requires small quantum spins, 
can be excluded from the beginning. 

Hence, we model the SWNT by a ID nuclear spin 
lattice of length L coupled through the hyperfine inter- 
action to a ID electron gas. The Hamiltonian resem- 
bles that of a Kondo lattice H — H e i + A^ i S, • I», 
where i runs over the ID lattice sites with positions 
r-j, lj = {If, If , If) is the effective nuclear spin of size 
I = IN±, Si = {Sf, Sf, Sf) is the electron spin operator 
at site i, and A = Aq/N± is the on-site hyperfine inter- 
action constant Aq weighted by the transverse electron 
mode. In contrast to the usual Kondo lattice model, H e i 
describes the interacting electrons and is defined in Eq. 
© below. 

The precise value of Aq in SWNTs is unknown. Es- 
timates in the literature [23[ provide values of Aq ~ 
1CP 7 — 1CP 6 eV, depending much on the curvature of 
the nanotube (higher values have been reported in (l7j 
though). This compares with the typical energy scales 
of the electrons, which can be quantified by the value 
Ep = vpkp/2 (we set H = 1 throughout this paper), 




FIG. 1: Illustration of the helical nuclear magnetism (indicated by the blue ribbon) of the single wall 13 C nanotube (SWNT), 
which emerges below a critical temperature through strong renormalization of the hyperfine coupling by electron correlations. 
The nuclear spins (red arrows) order ferromagnetically on a cross-section of the SWNT and rotate along the SWNT axis with a 
period n/kF in the spin xy-p\a,ne (chosen here arbitrarily orthogonal to the SWNT axis). Through the feedback of the nuclear 
magnetization the electric conductance of the SWNT is reduced by a factor of precisely 2. 



where k F /ir = n e i is the electron density in the system 
and v F (» 8 x 10 5 m/s in SWNTs [H, EE HI) is the 
typical velocity of electron excitations. Through the de- 
pendence on n e i, Ep can vary between the meV to eV 
range. 

Effective model. — Due to the small ratio A/ Ep, the 
energy and time scales related to the electrons and nu- 
clear spins decouple, and we can treat both subsystems 
separately. A Schrieffer- Wolff transformation of H allows 
us to obtain an effective Hamiltonian for the nuclear spins 
01, 
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where a = x^VjZ, and J™- = A 2 x?fa/2 is the ef- 
fective RKKY [25] interaction between nuclear spins. 
a is the lattice spacing and provides the short dis- 
tance cutoff of the continuum theory. The sum over 
q = mr/L for integer n runs over the first Brillouin 
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infinitesimal i] > 0) is the static electron spin suscep- 
tibility. We also have defined 1° = Yli e lri9 If and 
/ dr e~ lr iJ a {r). 
The effective electron Hamiltonian, on the other hand, 
includes the effect of the feedback of the nuclear field on 
the electrons. Since the spins I — IN± are large, we can 
choose H*f f = H e i + H 0v , with H 0v = J2i h » ' S* and 
hi = A(Ii) the nuclear Overhauser field. 

Interacting electrons as Luttinger liquid. — We use a 
bosonized Hamiltonian to describe the interacting elec- 
tron system of the armchair SWNT, which is naturally 
in the Lut ting er liquid state due to the linear electron 
dispersion |l4jl5| . The unit cell of a graphite sheet con- 
tains two carbon atoms, which results into a two-band 
description of the bosonized system. Since mixing be- 
tween the bands is essentially absent [lj, EH] we shall, 
however, focus on a single band only in order to avoid a 
heavy notation. The bosonized single-band Hamiltonian 
reads QUI, Hi 



v v K v {V0 v (r)f 
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where </> C;S are boson fields such that — V^csV^/tt ex- 
press charge and spin density fluctuations, respectively. 
9 C S are such that V0 CjS /7r are canonical conjugate to 4> c , s - 
v cs = v f/Kc,s are charge and spin wave velocities, and 
K c s are the dimensionless Luttinger liquid parameters. 
For SWNTs [3 El, K c « 0.2 . If the electron spin SU(2) 
symmetry is maintained, K s = 1, otherwise K s ^ 1. 

Without feedback from nuclear magnetic field. — Let 
us first assume that there is no feedback from the Over- 
hauser field on the electrons and set h.; = 0. The elec- 
tron system forms a Luttinger liquid, for which the zero 
temperature spin susceptibility has a singularity at mo- 
mentum q = ±2kp induced by backscattering processes 
[H, H3- At T > this singularity turns into a steep 
but finite minimum: The backscattering part of the spin 
operator Sf is expressed in the bosonization language by 
the operators (H O x SDW {n) cx e -^n e i^ c C0S (V26g, 
such that S x — [0% DW + Og DW ]/2 plus forward scatter- 
ing terms. Similar expressions [261 ] hold for S v and S z . 
We further assume that = J q is isotropic and in par- 
ticular K s — 1. The correlators between those operators 
can be evaluated in the standard way and we obtain (for 
q>0) 

J q (9, vf) w -C(g, v F ){k B Tf3- 2 \T(k)/T(k + 1 - g)\ 2 , 

(3) 

where g = (K c + K- 1 )/2, k = g/2-i\ T {q-2k F )/AT:, de- 
pending on the thermal length Xp — vp/ksT with k B the 
Boltzmann constant. Y is Euler's Gamma function and 
C(g,v F ) = A 2 asm(irg)Y 2 (l - g)(2na/vp) 2 a- 2 /8ir 2 v F - 
We have made the inessential assumption v c = v s = v F . 
Note that J q is independent of k F for a linear dispersion. 
A density dependence of J q requires a curvature of the 
electron dispersion, which partially restores Fermi liquid 
properties [28| , a scenario which we disregard for metallic 
SWNTs. A sketch of J q is shown in Fig. H 

At temperatures T < T * [defined in Eq. ([6]) below], 
\J2k F (T) | > k B T and the nuclear spins can - classically - 
minimize the RKKY energy by aligning in a spiral order 
I- = IN±[cos(2k F ri)e x + sm(2k F ri)e y ], where e X:V are 
orthonormal vectors defining the spin (x, y) plane. We 
shall henceforth assume that this order is established, 
and show that this assumption is self-consistent. Fluc- 
tuations reduce this maximal polarization, and in gen- 
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FIG. 2: Sketch of the RKKY interaction J q given by Eq. © . 



eral < IN±. The lowest lying excitations (to order 
1/IN±) in the nuclear spin system are magnons. Since 
Jij is long-ranged the energy cost of local defects, like 
kinks, scales with the system size and is very high. 

For a helimagnet, there exists a gapless magnon band 
with the dispersion @ w q = 2{IN 1 _)(J 2kp+q /N\ - 
J2k F /N\). Let rrii = (If) ■ l' i /(IN±) 2 measure the com- 
ponent of the average magnetization along IJ, normalized 
to < rrii < 1. Its Fourier component rri2k F acts as an 
order parameter for the spiral phase. Magnons decrease 
this order parameter and we have 0] 



m 2 k F (T) = 1 



E 
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where the sum represents the magnon occupation num- 
ber. In the continuum limit L — ► oo the integrand is 
divergent as 1/q 2 for q — ► (the q = mode is ab- 
sent because the system is not a ring), showing the ab- 
sence of true long range order in the ID system. Despite 
its appearance the divergence is not a consequence of 
the Mermin- Wagner theorem (2^, [3(3] , which forbids long 
range order in low-dimensional systems for sufficiently 
short ranged interactions. Since Jy is long ranged this 
theorem does not apply. 

The present situation, however, is very different in 
that the system has a finite length L ~ 2 /im imposed 
either through the natural length of the nanotube or 
through an external confining potential. At tempera- 
tures T < Tq we find that L < At, and so the cost of 
exciting the first magnon is already very high 0J q = n /L ~ 
2I\ J2k F (T)\/N±. We can define a temperature Tmo pro- 
viding the scale of the excitation of the first magnons 
by imposing u q /k B T » 2I\J 2kp {T)\/N ± k B T = 1. For 
T > Tmo we can then simplify Eq. ^ to 



m 2 k F (T) w 1 



1/IN± 
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where we have defined 

k B T* = [2I 2 C(g, v F )T 2 (g/2)T- 2 (l - g/2)) ^ . (6) 

For the SWNT this temperature satisfies the self- 
consistency condition HbTmo < k B T * <C vf/L. We use 
Tq as an estimate for the critical temperature. For a typ- 
ical SWNT Tq is very low. With the values given with 



Fig. [3] we obtain Tq ~ 10 /xK, too low for experimental 
detection. Yet this analysis completely neglects the feed- 
back of the magnetic field on the electron gas. This leads 
to a strong renormalization of T * . 

Feedback of nuclear magnetic field on electrons. — The 
ordering of the nuclear spins leads to a spatially oscillat- 
ing Overhauser field = that acts back on the 
electrons. We choose the electron spin axis such that 
S ■ e x — S x and S ■ e y — S v . The spatial oscillations of 

±2ik F Ti 



hi cx e 

tial oscillations of the Og'^ w 



in Hov perfectly cancel some of the spa- 
operators of the 5f ' y . Ne- 
glecting the remaining (irrelevant) oscillating terms we 
obtain H<jv ~ Si A Im2k F cos(v / 2~KV>+(ri)), where we 
have introduced 4>+ — (<j) c + 6 s )/yK with the normal- 
ization K = K c + 1/K S . The Hamiltonian becomes of 
the sine-Gordon type and Ho v is relevant in the sense of 
the renormalization group (RG): The 0+ field is pinned 
at a minimum of the cosine term of Ho v ■ The result is 
a density wave that combines charge and spin degrees of 
freedom. Fluctuations about the minimum are massive, 
with a mass associated to an energy scale A. At commen- 
surate electron filling Umklapp processes would become 
relevant too, and lead to fully gapped charge and spin 
sectors. For SWNTs, however, this would require high 
electron densities leading to Ep ~ 1.4 eV. This case is 
not considered here. 

Within a perturbative RG approach we find that 



(A Im 2 k F /E F ) 



V(2-ff) 



Vf I Ob. 



(7) 



This mass gap A is the first important consequence of the 
feedback. The second important consequence is the spon- 
taneous generation of anisotropy because the spin (x, y) 
plane is singled out by the Overhauser field. This is seen, 
for instance, in the spin susceptibilities x aa ■ Those can 
be calculated in the same way as before (see Appendix A) 
if we notice that the massive ip + field does not contribute 
to the long-wavelength asymptotics. The finite temper- 
ature expressions for the x aa are otherwise identical to 
the case without feedback, and the RKKY couplings Jg 
can be obtained from Eq. ^ upon the following modi- 
fications: For x xx an d X vv the exponent g is replaced by 
g' = 2K C /K S K and the amplitude is reduced by a factor 
2 because a term depending on <f> + only drops out. For 
X zz the exponent becomes g" — (K c /K s + K C K S )/2K 
while the amplitude remains unchanged, vp is replaced 
by v- = (v c /K c + v s K s )/K. This leads to 



J*>« = J q (g',V-)/2, JJ 



(8) 



For K c = 0.2 and K s = 1 we have to compare g = 0.6 
with the strongly renormalized g' = 0.33 and g" = 0.17. 

Let us finally note that correlators between <f> + , 9 + can 
only be neglected as long as k B T < A, i.e. A^ 1 < 
with £ = vp/A the correlation length. In Eq. © below 
we define a critical temperature T* similarly to Tq before. 
For T < T*, m 2kF w 1 (see Fig. [3]), and we find that 
A ^> k B T. At T — > T* , however, m 2 k F vanishes and 
so does A. The order in electron and nuclear systems, 
therefore, vanishes simultaneously. 
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Consequences for magnetization and transport. — The 
helical order still minimizes the energy and there remains 
a gapless magnon band H, uj q = 21{J' 2kF + q 
where J' q = J q 

magnetization remains otherwise unchanged. Replacing 
J q by J' q in Eq. ([6]) leads to the renormalized critical 
temperature T*, 



jy. The previous discussion of the 



k B T* = [/ 2 C(.g',t-_)r 2 ( 5 '/2)r- 2 (l -g'/2)] 



(9) 



The notable difference is the modified exponent. For the 
parameters displayed with Fig. [31 we obtain the change 
from 1/(3 - 2g) = 0.625 to 1/(3 - 2g') w 0.43. Quite re- 
markably this considerably boosts the value of the char- 
acteristic temperature from Tq ~ 10 to T* ~ 1 
mK. Note that T* -C vp/Lks is still satisfied. Fig. 
[3] (a, solid line) shows the result of the feedback. In 
Fig. [3] (b) we also show the dependence of T* on Aq, 

,2/(3-2,,') 
A 



T 



= Aq- S6 



The order furthermore modifies the transport proper- 
ties of the system. With the opening of the mass gap in 
the (p+ channel, half of the conducting modes are blocked 
and the conductance decreases by the universal factor of 
2. As an illustration we consider a SWNT connected to 
metallic leads. The conductance is given by [HI, HH, [33| 
G = 4e 2 //i, where e is the electron charge, h the Planck 
constant, and where 4 is the number of conducting chan- 
nels (2 spin projections and 2 bands). The pinning of the 
(j>+ field (in each band) blocks 2 conductance channels 
and so reduces the conductance precisely by the factor 
2 (see Appendix B for details). Such a reduction is a 
direct consequence of the nuclear spin ordering and the 
Luttinger liquid physics of the electrons, and should be 
detectable experimentally in standard transport setups. 

As a conclusion, we emphasize that the physics de- 
scribed here is quite general and is also relevant for other 
ID systems of the Kondo-lattice type. 

This work was supported by the Swiss NSF and NCCR 
Nanoscience. 



APPENDIX A: SUSCEPTIBILITIES IN THE 
PARTIALLY GAPPED SYSTEM 

The static susceptibilities x a ( a ) = X Q (<Z: W = 0) are 
the Fourier transforms of the following retarded electron 
spin response functions (we set h = 1 throughout this 
section) 



X a (r,t) = -te(t)({S a (r,t), S"(0,0)]>, 



(Al) 



where S a (r, t) are the electron spin operators at position 
r and time t, and Q(t) is the step function. 

In the bosonization language, the spin operators are 
expressed through the sum of forward scattering and 
backscattering operators as [26[ S a — {S a )f orw + 
[Osdw + i^SDwV]/^- The forward scattering opera- 
tors lead to a contribution to x° that is regular in the 



momentum q, while the backscattering parts are singu- 
lar or strongly peaked at momentum q = 2k p. We shall 
therefore neglect the forward scattering contribution and 
focus only on the backscattering operators (not writing 
Klein factors) [H 



O 



SDW 



W SDW 



o 



SDW 



p -2ik F r 

2na 

„-2ik F r 

i- e 1 ^ 

2na 

„-2ik F r 

i- 



2na 



is/28 3 



iV20 s 
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(A2) 
(A3) 
(A4) 



The average in Eq. (|Alj) is evaluated with respect to the 
following Hamiltonian, depending on the gapped (f> + and 
ungapped <fi- fields, 



H 



[(VMr)) 2 + (VMr)^ 
(V0 + (r)) 2 + (V0 + (r)) 2 +r 2 <Pl(r)]}, 



(A5) 



where £ = v + /A is the correlation length associated to 
the gap A, v + = (v c K c + v s /K s )/K and w_ = (v c /K s + 
v s K c )/K for K = K c + l/K s . 

The boson fields <f>± , 9± are related to the usual charge 
and spin boson fields through the transformation 



(A6) 
(A7) 
(A8) 
(A9) 




Let us set r — (r,t), x a { r ) — — *2(a7r) 2 x a (?", t), and as- 
sume that t > 0. We can then write, for instance for 

x x , 

x x (f) = cos(2fc F r) \([ e ^<t> + (r) , e -iv^0 + (o) j ) 

(A10) 

with K' = AK C /K S K and K" = (K c - K^f/K, and 
where we have used the invariance of the Hamiltonian 
under a simultaneous sign change of all the boson fields. 
In the Gaussian theory (| A5|) the correlators in (|A10[) 
are fully expressed through boson correlators of the form 
(0±(^)0±(O)) which, in (q,tu) space, are given by [26j 



((f)*_(q,uj)4)-(q,uj)) = 

for the massless fields and 
{(f>\{q,uj)(j) + {q,uj)) 



(u> ± irf) 2 — v 2 q 2 ' 



(uj ± ir\) 2 — v\q 2 — A 2 



(All) 



(A12) 
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FIG. 3: a: Magnetization m,2k F (T) [Eq. Dashed line: without feedback. Solid line: with feedback. Parameters for the 

curves are Q El HI HI E F = 0.1 eV, A = 10" 7 eV, v F = 8 x 10 5 m/s, a = 2.46 A, K s =l,K c = 0.2 (leading to g = 0.6, 
</ = 0.33), and L — 2 /xm. The vertical lines mark the temperatures written next to them, b: Characteristic temperature T* 

[Eq. ([9]l] as a function of the hyperfine constant Ao. The curve follows a power law T* oc A /{i 9 ' 
to the self-consistency limit T* w vf/LIcb =3K. 



aO.86 
A , 



and is plotted up 



for the massive fields. The sign the infinitesimal shift ±irj 
is dictated by the time order of the operators. Eq. (|A11[) 
is singular at to — > 0, q — > 0, and the proper treatment of 
this singular behavior leads to the power law behavior of 
the susceptibilities (at the shifted q — > ±2kp ) character- 
istic for a Luttinger liquid theory. 

On the other hand, Eq. (|A12[) is regular at w — ► 0, g — > 
0, and so does not contribute to the power law divergence 
at q — > ±2kp. In fact, let us expand Eq. (|A10|) in powers 
of c/>+ . The lowest nonzero term in the c/>+ is 



cos(2fci?r) 
+ 2Jf"([e 



2if([0 + (f),0 + (O)]) 

J™*- & 4>+ (f) , e" 4 ^ 7 *- (°) 0+ (0)] ) 



(A13) 



The Fourier transform of the first term, depending on 
4>+ only, at w — > 0, q — > ±2kp tends to a constant oc 
1/A 2 , and so contributes only insignificantly to the static 
susceptibility x x (l, UJ = 0) at q « 2fcp. The second term 
involves a sum of products of the type 

cos(2fc F r)e K ' (0 - (f)< ^ (0)> (<?!>+ (f)0+ (0)) . (A14) 

Since <f>- is massless, the exponential factor evaluates to 
a power law. With Eq. gig) we can then write the 
Fourier transform of Eq. (|A14[) as a convolution of the 
form 



dq'cLu' 



1 



v 2 _q' 2 



l-K' 



where q± — q ± 2kp. For cj — > and (j± 



q') 2 - A 2 ' 
(A15) 
we see that 



the u>' integral is dominated by the poles at iyq" 



q' 2 + A 2 

and the weak singularities at ±q'. The contribution at 
the poles is more singular, and if we focus on the pole at 



v\q' 2 



A 2 , we obtain 



1 



A 2(l-K') 



dq 



v\q n 



A 2 



(A16) 



The remaining integral leads to an arcsinh, which has 
an ultraviolet divergence that has to be cut off at l/a. 
More importantly, however, the result has no infrared 
divergence, meaning that this expression remains regular 
at q± — > 0. 

The latter results allow us to conclude that the Fourier 
transform of Eq. (|A13[) is regular at q — > ±2kp. Since 
the theory (|A5|) is Gaussian, higher order correlators are 
products of the latter results and so remain regular. We 
have therefore shown that the singular behavior of the 
susceptibility is fully controlled by the (/>+ independent 
term in the expansion of the e ^ + , allowing us to use 
the approximation 



X x (r) 



<2k F r)( 



]>> 

(A17) 



which is of precisely the same form as the susceptibility of 
a regular Luttinger liquid. The difference in the present 
case is the modified exponent K' and the fact that the 
first term in Eq. (|A10|) . depending on <ft + only, drops 
out. The amplitude of the resulting susceptibility X x {q) 
is reduced by a factor of 2. With g' = K' /2 this leads to 
the zero temperature susceptibility [26| 



X x (q) 



1 sin(7r#') ,, 



2 Av. 



r 2 (i-<?') 



± 



a(q±2k F ) 



2-2g' 



(A18) 

The susceptibility x v (q) involves the same combinations 
of the </)+ and cf>- fields and is identical to x x (q)- The sus- 
ceptibility x z (q) is expressed through the combinations 

± 4>, = -fe. (K c 4> + - x ff(f>- T&+T VKJT S 9-). In 



l 

contrast to x x an d x V i the massless fields do not cancel 
out for one of the ± signs, and so the previous reduc- 
tion of the amplitude by the factor 2 does not occur. 
Again we can neglect the contributions from the massive 
<f> + and 0+ fields. The resulting exponent of the power 
law in x z {q) then depends only on the combination of 
the prefactors of the (j>- and 6L fields, and is given by 
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(K c /K s + K C K S )/2K so that where g r (r, r'; us) is the Fourier transform of the retarded 

boson Green's function 



a(q±2k F ) 



2-2g" 



g r (r,r ] t) = -ie(t)(lq> c (r,t),Mr',0)}) 
( Al9 ) =L + (ry)g r + (ry;t) + L.(ry)g r _(ry;t), (B2) 



The extension to temperatures T > is straightforward 
(2^ | and the result is given by Eqs. © and ©. The 
main effect is that the singularity at q = ±2fcp turns with 
into a finite minimum. 

As temperature rises, the depth of this minimum de- g±(r,r;t) = — i0(t)([0±(r, t) , <f>±(r',0)]) (B3) 

creases. The neglected <fi+ correlators become important 
when the temperature becomes comparable to A. For 
temperatures below T* , the helical magnetization m2k F anc ^ 
is close to 1 and so A w {Al j 'Epfl^^vp /a. For the 

values (see Fig. \5§ A/E F = lO^ 6 ,^ - 8 x 10 5 m/s, L , f) = K c {r)K c {r') 

a = 2.46 A and 5 = Jf/2 = 0.6, we obtain A w 6.8 x +K r ^ r ! ^/K{r)K{r') ' 

eV, i.e. a corresponding temperature of 0.7 K. This is 



much higher than the critical temperature T* and con- j , i\ 

K c {r)K c {r>) ,_ c . 
firms the validity of the approximations above. ' y K s {r)K s {r')K{r)K{r') 



APPENDIX B: REDUCTION OF THE 
ELECTRICAL CONDUCTANCE 

We illustrate here the reduction of the conductance of 
the nanotube by a universal factor 2 with the calculation 
of the electrical DC conductance G of a Luttinger liquid 
that is connected to metallic leads. It was shown in [3ll . 

I32I I33I ] that in such a system the conductance is given by 

2 

G = j-n, where n is the number of conducting channels 
in the Luttinger liquid. The Luttinger liquid properties 
appear in this formula only through the number n of 
channels, while the prefactor e 2 /h is entirely determined 
by the properties of the leads. 

A metallic single wall carbon nanotube has n = 4, 
which is composed of a factor 2 arising from the two spin 
directions, and a factor 2 because the unit cell has two 
carbon atoms leading to two electron bands. With the 
feedback from the ordered nuclear magnetic field the </>+ 
field is pinned and we expect that it no longer contributes 
to the electrical conduction. Since there is a <p + in each 
band, the number of conducting channels then reduces 
to n = 2. 

In the following we shall prove this intuitive result: As 
in [3J HH, [33| we model the leads by one-dimensional 
Fermi liquids with K c = K s = 1. The validity of the 
quadratic Hamiltonian (|A5[) can then be extended into 
the leads by introducing a spatial dependence on v± and 
£ (or A) such that v±{r) = vf and £(r) = A(r) = when 
r lies in the leads. 

In the bosonization formulation the DC conductance 
can be calculated through the u) — > limit of the nonlocal 
conductivity [26[ 

4e 2 

a{r,r;u) = -i—{uj + ir])g r {r,r';Lu), (Bl) 



As found in [3l|, [32|, [33j the conductance depends only on 
the properties of the leads, where L±{r,r') = 1/2. If the 
system were entirely gapless, both and 0+ channels 
would therefore contribute equally to the total conduc- 
tance, with 2e 2 /h from each channel (the latter factor 2 
being due to the two bands). 

The gap A affects only the <fr+ fields, and so the con- 
tribution 2e 2 /h from the two 4>- channels remains un- 
changed. The contribution from the gapped cf>+ fields, 
however, drops to zero: If we follow the method of Ref. 
[3l| the inclusion of the finite gap A(r) in the calculation 
of g r + is straightforward. The result is 



g r + {r,r';uj) = 



1 



3 _| r - r '|VA3 



2VA 2 



(B6) 



not writing further regular terms expressing reflections 
of propagating waves on the ends of the nanotube, and 
for r and r' both lying in the nanotube. In the DC limit 
ui — > we then see that the corresponding conductance 



{u + irj)g r + {r, r';u) -> 0, 



(B7) 



while ungapped fields have here the nonzero limit i/2. 
Hence the DC conductance is entirely determined by the 
two 4>- channels from the two bands of the nanotube and 
is given by 



(B8) 
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